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The coherent potential approximation, CPA, is a useful tool to treat systems with 
disorder. Cluster theories have been proposed to go beyond the translation invariant single- 
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site CPA approximation and include some short range correlations. In this framework 
one can also treat simultaneously diagonal disorder (in the site-diagonal elements of the 
t-h , Hamiltonian) and non-diagonal disorder (in the bond energies) . It proves difficult to obtain 
reasonable results, free of non-analyticities, for lattices of dimension higher than one [D > 1). 
' We show electronic structure results obtained for a Hubbard model, treated in mean field 
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approximation, on a square lattice and a simple cubic lattice, with the simultaneous inclusion 



of diagonal and non-diagonal disorder. We compare the results obtained using three different 
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methods to treat the problem: a self-consistent 2-site cluster CPA method, the Blackman- 



Esterling-Berk single-site like extension of the CPA and a nested CPA approach. 
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Many-site or cluster extensions to the single site CPA introduced by Soven in 1967 jjj have 
been considered (see Refs.[2], e.g.) to overcome its limitations. Few results are available for 
D > 1, while many approaches have limited applicability due to spurious non-analyticities 
in the Green's function ||. Analyticity of the theory was proved in a few cases ||-^|. CPA 
becomes exact in the infinite-D limit , and interesting connections can be made with the 
recently developed methods for strongly correlated systems ||. 

We now briefly present our electronic structure calculations (a more complete presentation 
will be given elsewhere |))). The Hamiltonian considered here is: 

H = E E c4°j + E E * a t a i - 1 E c t c j 

jec jga <t,j>ec 

-t E 4 a s-* E • // <- (!) 

<i,j>£A <i,j>,ieA,jeC 

Cj(Oj) creates an electron in the C(A)-species tight-binding (TB) orbital centered at site 
j. The two components of the alloy, a = C, A respectively, are described by a single 
band Hubbard model JlC| ] each, with site energies e a , local correlation U a , and hopping 
t. A mean field treatment for the paramagnetic phase of the model maps the problem 
to a TB Hamiltonian with renormalized site energies, here: Ec = £c + Uc <T 2 C> > an d 
Ea = £a + Ua <U 2 > ■ describes hopping between sites occupied by different species. 

To treat this model we first employed a self-consistent 2-site cluster CPA approach (de- 
noted hereafter 2-CPA) [pd|Js|] . A uniform effective medium is defined, through a 2x2 complex 
self-energy matrix S to be determined self-consistently by the generalized CPA-condition 
for the cluster. Concretely, if Ho denotes the tight-binding Hamiltonian of a D-dimcnsional 
lattice with hopping £, we define: 

Heff = H o + £ , Sjj = {<Jd - a n ) <5i,j + cr n - (2) 

Isolating a cluster composed of two nearest neighbour sites immersed in the above effective 
medium, one has the reduced cluster Hamiltonian: 

H(i,2) = H eff (cj) + A#H, (3) 

where AH (lj) is the 2x2 matrix describing the difference between the projection of the real 
Hamiltonian H onto the two sites of the cluster, and the effective medium. We impose that 
the average Green function should coincide with the effective medium Green function, G e //. 
Using Dyson's equation, this "2-CPA" condition can be written: 

< [I- AHGeff^AH >(i, 2) = 0. (4) 

The second method we employed is a variation of the 2-CPA, the nested CPA approach 
(NCPA). Basically, one can rewrite Eq. Q in a convenient form such that one can split 
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the problem of finding the self-energy E, in many (five) simpler CPA alloy calculations, 
which are "nested" as three consecutive levels of CPA are required. The average over the 
occupations of the 2 sites of the cluster implied by Eq. Q has the form || : 

Pec Xcc ~t~ Paa Xaa "t - Pac Xac Pea Xca 0> (5) 

where c is the concentration of the A-species, p is a parameter included to analyze phase 
segregation, and: P cc = (1-c 2 ) + c(l-c)p, P aa = c 2 + c(l-c)p, P ac = P ca = c(l-c)(l-p). 
Introducing Xt = (Xac + Xca)/2, it can be rearranged in a form suggestive of the sequence 
of alloy problems to be solved in the NCPA || : 

= P [(1 - C) X cc + C X aa] + (1 -P) {(1 " C) 
[(l-c)Xcc + C X t] + c[c X aa + {1 - c) X t}} ■ (6) 

Finally, the third method we employed was the single-site extension of the CPA by Black- 
man, Esterling and Berk (BEB) f| to include non-diagonal disorder. 

In Fig. 1 we exhibit the results we obtained with the three methods described on a 
square lattice at c=0.5, for which we used the exact bare local and nearest neighbour's 
Green functions in terms of the complete elliptic integral (T^j, continued analytically onto 
the physical sheet of the Riemann surface || . The corrections to the single site approach of 
BEB obtained with the 2-CPA and the NCPA are not qualitatively very important, except at 
the band center. In Fig. 2 we show results for a three dimensional lattice: here we adopted 
the semiclliptic approximation to the bare Green functions of the simple cubic lattice. 
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Figure Captions: 

Fig.l - Square Lattice: Density of states (DOS) vs. energy, approximations employed 
labeled inside plot. Energies in units of the hopping t. Open symbol: C-species local 
density of states in 2-CPA approx. Parameters: Ec = 0, £U=2.5, t'/t=l.l, c=0.5, p=0. 

Fig.2 - Simple Cubic Lattice: DOS vs. energy. E c = 0, E A =2, t'/t=1.2, c=0.25, p=0. 
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Fig. 1 - C.I.Ventura - SCES Ref.# 373 
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